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Abstract: This paper explores whether quantum field theory allows the events of emission 
and absorption of a single particle to be separated by a space-like interval without violating 
Lorentz symmetries and causality. Although the answer is indeed affirmative, tradition- 
ally such effects have been considered negligible. We show that for sufficiently light mass 
eigenstates such processes can become significant over macroscopic length scales. A critical 
review of the historical literature reveals various shortcomings of the standard methods; 
specifically, one finds that they are restricted to states for which the expectation value of 
momentum vanishes. Furthermore, the results obtained here correct Feynman's analysis 
of this subject. A formalism is thus developed that allows the description of states with 
non-zero momentum, which is then applied to the OPERA and ICARUS neutrino-speed 
experiments. For OPERA we choose a mass in the nano electron-volt range and find 
that although our formalism predicts a non-zero detection probability for an early arrival 
time of 60 ns, the predicted event distribution is maximal on the light-cone. Consequently, 
our prediction does not reproduce the peak at 60 ns reported by the OPERA collabora- 
tion. Turning to the ICARUS experiment, we note that while the collaboration reported 
an average time of fiight that is consistent with the speed of light, the event data with 
its associated uncertainties nevertheless indicates that some of the detection events are 
separated from their corresponding emission events by a space-like interval. For a micro 
electron-volt mass range, this is in agreement with the here reported formalism. We thus 
raise the possibility of employing high-precision neutrino-speed experiments to determine 
the absolute masses of neutrino mass eigenstates. Finally, it is argued that the predicted 
space-like amplitudes are consistent with the emperical data on SN1987a. 
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1 Introduction 

This communication is a preliminary exploration of the possibility that massive particles 
may affect events at space-like separations. We begin with the observation that in the 
absence of quantum mechanical considerations Lorentz, symmetries unambiguously divide 
space-time into causally connected and causally disjoint regions. The boundary is occu- 
pied by massless particles. In classical theories, no massive particle can trigger a detector 
placed at a space-like separation from its source. In quantum theories, however, the non- 
commutativity of the position and momentum operators allows a classically-forbidden trig- 
gering of detectors. This possibility, as far as we know, first arose in the work of Pauli [1]. 
He was troubled that this could lead to a violation of the causal structure of space-time 
[1-3]. It was not until the work of Stiickelberg and Feynman [4-6] that the path for a 
consistent resolution was opened: for such processes the time ordering of events is not 
strictly respected, yet causality is preserved by the introduction of antiparticles. For the 
frames in which temporal ordering is reversed, what one would have seen as the absorp- 
tion of a particle is instead observed as the emission of an antiparticle and vice versa. At 
the textbook level, this resolution is presented in the monographs by Feynman [3, chapter 
18] and Weinberg [7, section 2.13]. The understanding commonly conveyed in the mod- 
ern literature is that the effects are too small to have directly observable consequences at 
macroscopic scales, see for example [8]. 
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This conclusion is somewhat misleading because the scale inextricably depends on the 
mass of the particle. For example, Weinberg's analysis shows that the scale at which such 
an effect becomes significant for a proton is roughly 2 x 10~^^ cm. Nevertheless, it seems 
to have been overlooked that if one goes down in mass by eighteen orders of magnitude 
to the nano electron-volt (neV) scale, then even in the Weinberg-like back-of-the-envelope 
calculation the range of the said effect immediately stretches to 2 x 10^ cm, that is, to 
two hundred meters. Given that the present data on neutrino oscillations places no lower 
bound on the absolute masses of neutrino mass eigenstates [9] , this observation raises the 
possibility of novel neutrino phenomena. Indeed, our analysis suggests that neutrino-speed 
experiments may be employed to experimentally measure the absolute masses of neutrino 
mass eigenstates provided one of them is sufficiently light. 

Even though this observation was inspired by the OPERA neutrino-speed experiment 
which has been recently called into question, the here-presented analysis remains valid, and 
may be of importance for a wider class of neutrino-speed experiments. One such example 
specifically discussed in this communication is the ICARUS experiment. 



2 Amplitudes for space-like separations and causality 
2.1 Calculation of the amplitudes 

To give a concrete form to the above discussion we must calculate the amplitudes for the 
classically-forbidden triggering of detectors mentioned in the opening paragraph of this 
communication. Rather than following the simplest example of a scalar field, we calculate 
these amplitudes in the context of a spin one-half Dirac field. This approach allows for the 
discussion of a number of conceptual questions that would otherwise remain unexposed. 
We thus consider the field [10] 

^K^) = E / 7^ a/^ h ^) e"'""" « (P' ^) + (P' ^) ' (2.1) 
{IttY'^ V p l J 

where the index i, £ G {1,2,3,4}, labels the four components of the well known Dirac 
field ^'(x).^ Viewing equation (2.1) as a linear combination of creation and annihilation 
operators, one naturally arrives at the question: what does the operation of ^*{x) on |vac) 
yield?^ In the analogous case of a scalar field, Hatfield interprets the resulting object as 
a state localised at x [11]. In the context of equation (2.1), such an object is in fact a 
four- dimensional column with components 

\x)i := *K^)|vac). (2.2) 

The dual to \x)i is 

^{x'\:=Y,{^sic\^e'{x'he'e, (2-3) 



^The notational details are given in appendix A. 

^The reader is referred to the remarks following equation (A. 8) for the definition of *. We consider the 
operation of '^*{x) on |vac) rather than that of '^{x) on |vac) purely for convenience. 
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where is introduced as a metric in order to yield a Lorentz invariant inner product. 
Therefore, the unnormalised amphtude for a particle created at x to trigger a detector 
located at x' reads^ 



A{x x') := i{x'\x)i 
I 

= Yl (vac| (x')7° (x) I vac) 



= (vac|^''^(a;')7°^'*(2;)|vac). (2.4) 

The separation between x and x' may be time-like, light-like, or space-like. The last of 
these separations shall yield the classically-forbidden triggering amplitudes. Evaluating the 
right hand side of equation (2.4) using equation (2.1), one obtains 

_ f d^p 2m 



A similar evaluation of the amplitude for an antiparticle yields 

A{x -f x') = (vac|^'"f(x')7°^'(2;)|vac) 

= -A{x^x). (2.6) 

These amplitudes, modulo the minus sign in the above equation (dictated by the fermionic 
nature of the Dirac field), coincide with the known result for a massive scalar field. ^ 

The evaluation of integrals such as equation (2.5) is far from trivial, and was first 
undertaken by Dirac [12]. Following his approach, one can rewrite equation (2.5) in the 
form 

A{x^x') = -i^^U{r,t), (2.7) 
iTir or 



where^ 



1 f °° 

U{r,t) = — g-im[r-smh(vp)+tcosh(^)] ^2.^ 

TTi ./-oo 



Here ip is the rapidity parameter 

IpI = msinh(93) and = m cosh{(p) , (2.9) 

and 

r:=\x'-x\ and t := \x'° - x^\. (2.10) 

To proceed further, we employ an integral representation of the Hankel function of the 
first kind that can be identified with integral (2.8). By considering a Lommel's expansion 



^This interpretation appears implicitly in Hatfield's analysis of the problem [11]. 

*Tlie apparent spin-independence of A{x — > x') and A{x — >■ x'), apart from the indicated sign, is most 
likely due to the absence of any (spin-dependent) coupling to external fields. 

^The here defined U{r,t) differs by a sign in the exponent from the expression quoted by Dirac [12]; this 
is consistent, because Dirac was integrating a slight variation of equation (2.5). 
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[13, page 140] of the Hankel function of the first kind and equating it to a suitable series 
expansion, along with an appropriate change of variables, one obtains (in the case where 
the order v of the Hankel function of the first kind is set to = 0) 

1 /"OO , , 

_ / g-im[icosh(^)+irsinh(v)] ^ fjW 1 ^2 _ , (2.11) 

Mathematically there also exists an alternate solution to the aforementioned series ex- 
pansion in terms of the Hankel function of the second kind, however, this solution proves to 
be unphysical at space-like separations and equivalent up to a sign at time-like separations; 
consequently it is not considered. Substitution of equation (2.11) into equation (2.7) yields 

2 

A{x x') = -— , ^ (m^¥^\ . (2.12) 

27r _ \ / 

This expression is valid for all space-time separations. If x and x' are separated by a 
purely space-like interval, amplitude (2.12) may be rewritten in terms of a modified Bessel 
function of the second kind of order u = 1 

A(x x') = -A(x x') = ^—II^=Ki ( m^/r'^ - A . (2.13) 

The result (2.12) is consistent with that of Pauli [1], and for an analogous integral in a 
different context, also with that of Dirac [12]. Feynman, however, obtains a result in terms 
of a Hankel function of the second kind and order v = 1 [3, 6]. This result proves to 
be divergent for an increasing parameter of the space- like interval [13]. A second error 
occurs when he gives an incorrect asymptotic expansion for Hankel functions of the second 
kind [13, page 198]. For a space-like interval, the resulting incorrect expansion matches 
the asymptotic behaviour of equation (2.12); however, inside the light-cone it yields an 
oscillatory solution that contradicts equation (2.12). 

Of course the true genius in these historical investigations lies in the prediction of 
antiparticles. We take a brief moment to revisit this argument. First note that since 
^2 _ ^2 Lorentz invariant, the amplitudes A{x — >• x') and A{x — >• x') also exhibit the 
same invariance. Second, consider two space-like separated events x and x' . Since the 
amplitudes A{x — t- x') and A{x — t- x') are non-zero there exist two classes of observers, one 
for which x'^ > x^ and the other for which x^ > x'". The temporal ordering of events is 
thus no longer preserved. Causality is then ensured by the following interpretation: if the 
first set of observers sees a fermion created at x and subsequently absorbed at x', then the 
second set of observers sees an antifermion created at x' and absorbed at x, and vice versa. 
As far as we can infer, the credit for this insight goes to Pauli, Stiickelberg, and Feynman. 
The extension of these arguments to bosons runs along similar lines. 

2.2 The massless limit 

In this section we explore the behaviour of the amplitude (2.13) in the massless limit. As has 
been shown by Weinberg, the massless limit of the Dirac field is well defined and coincides 
with the field operator obtained by the demand that any Hamiltonian density constructed 



-4- 



from the field operator must transform appropriately under Lorentz transformations [14]. 
Consequently, we are justified in directly taking the m — )■ limit of the amplitude (2.13). 
For a finite space-time interval — the argument of the modified Bessel function of 
the second kind will be small as m approaches infinitesimal values, hence we can consider 
the limiting form 

KAz)\^^,-^l{iy-iy-(l:z) , (2.14) 



for z € C and Re(i^) > [15]. It thus follows that the amplitude for a classically forbidden 
triggering vanishes in the m — )• limit 

lim A{x x') = lim -^r—r ^ = 0. (2.15) 

m—>-0 m— s>o vr (r — t ) 

Furthermore, for a fixed space-time interval, the amplitude (2.13) vanishes for both large 
as well as small values of m with a maximal probability dictated by the order of the space- 
time interval. Figure 1 depicts the mass dependence of the unnormalised^ probability 
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Figure 1. The plot depicts the mass dependence of the probability V{x x') for a fixed space-time 
interval of ^/r^ — t'^ = 1 neV~^ = 1.24 km. It highlights the existence of a maximal probability with 
respect to a variation in mass and, consequently, the smoothness of the massless limit. 



V{x —7- x') := J^ix — >• x'')A{x — )■ x') for a space-time interval of V?"^ — = IneV ^. It 

^The probability is unnormalisable due to the highly localised nature of the state that results from the 
application of a field operator on the vacuum (see section 2.3 below). The remedy is, of course, well known. 
One must introduce a state with a finite width in configuration space. 



-5- 



is possible to further qualify this result: repeating the above derivation of equation (2.13) 
using a massless Dirac field leads to the same result as that obtained in equation (2.15). 

2.3 Interpretation of the amplitudes 

In order to provide further physical interpretation for the above amplitudes a closer look 
at |x) is required. From equation (2.2) this state in component form reads 



It is thus apparent that for \x) the expectation value of the momentum operator vanishes 
with infinite variance. This agrees with the expectation for a completely localised state 
in configuration space. Consequently, the amplitudes A{x — )• x') and A{x — )• x') are only 
applicable to particles for which the expectation value of momentum vanishes. 

The application of these amplitudes to neutrino-speed experiments [16-18] requires 
additional effort.^ Before we embark on this task in section 3, we note that a graphical 
analysis of the probability V{x — )■ x') reveals that a particle has the highest probability of 
triggering a detector located at a light-like separation from its source. Furthermore, one 
finds that there is an asymmetry between the probability of a particle triggering a detector 
at a space-like separation and a time-like separation. 

3 Neutrino-speed experiments 

We will now explore a context in which the above amplitudes have an experimentally 
observable effect; in particular, we consider neutrino-speed experiments. 

While neutrinos produced in electroweak processes are a linear superposition of mass 
eigenstates, we make the simplifying assumption that the entire process corresponds to 
the emission and detection of a single mass eigenstate. This assumption will capture the 
essence of our argument. Thus, consider a single ultra-relativistic mass eigenstate produced 
at a time rj = and directed towards a detector placed at a spatial distance L. In the 
absence of the discussed amplitudes, the canonical expectation is that it may be detected 
no sooner than tc = L. In the spirit of this preliminary investigation, at any given instant r/ 
the emitted mass eigenstate may be instantaneously considered at rest (this allows the use 
of the amplitudes obtained above). Its classically-expected spatial location is r/. One can 
thus ask: what is the probability amplitude that the mass eigenstate under consideration 
triggers the detector at a time to < tc? That is, what is the probability amplitude, A{T), 
that the mass eigenstate triggers the detector at a time T := tc — to earlier than the 

^The work of Morris [19] considers the evolution of a wave-packet to describe a moving particle in this 
context. Such a description is, however, not valid in a relativistic context; see, for example, H. and T. 
Padmanabhan [20, section V. A.] for a detailed discussion on the evolution of a relativistic wave-packet. 




(2.16) 




(2.17) 
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canonically expected time? Within the defined setting, an apphcation of equation (2.13) 
leads to the result 



^(r)= [W^^K.m], (3.1) 

Jo ^iv) 



where ^(ry) := mY^(ir^^7/)2'^^'(t^^'^^^^p. 

In order to circumvent the unnormalisable character of P{T) we define the ratio 

(r^ - Zill - -^*(r)-4(r) 
^''^ ' p{r,)~ A*{r,)A{r,) ^'^ ' 

where is a conveniently chosen experiment-specific reference point. 

3.1 OPERA neutrino-speed result 

For comparison with the OPERA results, we choose Feo •= 60 ns as the reference point, and 
use equation (3.2) to calculate a6o(r). For a fixed baseline L = 730km, as is applicable to 
the OPERA experiment, a detailed graphical analysis shows that as one lowers the mass in 
expression (3.2) the variation of ago with F saturates for m < 10~^ neV. This mass range is 
a natural choice in the context of the OPERA experiment because it maximises the space- 
like amplitudes in expression (3.1). A graphical representation of a6o(r) for m = 10~^ neV 
is provided in figure 2, which reveals that although there is a non-vanishing probability 
of detecting a mass eigenstate at F ~ 60 ns the predicted event distribution is larger for 
smaller values of F.^ The most expedient path to compare this event distribution with the 
OPERA results is to draw attention to the data pertaining to the short-bunch wide-spacing 
beam [16, figure 18]. For the 20 events of this type recorded, the collaboration reported 

a. a mean value of F = (62.1 it 3.7) ns, 

b. that all the 20 events occur roughly in the range 40 ns < F < 90 ns, 

c. that there are no events in the range for F < 40 ns. 

This does not match our theoretically predicted event distribution. 

Since the original announcement, the OPERA collaboration has unearthed two sources 
of errors [24] . At the time of writing it is not clear how the proverbial dust will settle [25-39] . 
Nevertheless, provided that there exists a sufficiently light mass eigenstate, our analysis 
suggests that the space-like correlations under discussion will be detectable. 

3.2 ICARUS neutrino-speed result 

We now turn to the neutrino-speed experiment performed using the ICARUS detector. 
Although the collaboration reported an average time of flight consistent with the speed 
of light, the individual events detailed in reference [18, table 1] nevertheless indicate that 
some of the detection events are separated from their corresponding emission events by a 



Integral (3.1) was solved numerically by employing the mpmath arbitrary-precision floating-point arith- 
metic library [21]. Figures throughout this communication were produced using the matplotlib plotting tools 
[22] in conjunction with the Numpy and Scipy scientiflc computation libraries [23]. 
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Figure 2. A graphical representation of the probabihty to detect a single mass eigenstate outside 
the light-cone. The ordinate ago corresponds to the ratio of the probability at T and the probability 
at Fgo := 60 ns. 



space-like interval. It is the distribution of such events that we wish to compare to the 
here-predicted event distribution. 

Guided by the ICARUS data, we choose Fio := 10 ns and compute the distribution 
aio(r). To gain a qualitative insight, the results obtained for a representative mass of 
m = 1 iieV are depicted in figure 3. Although this choice of mass roughly reproduces the 
ICARUS distribution [18, figure 3], it is not possible to perform a statistically meaningful 
fit given that the currently available data consists of only seven events with a subset of 
four events corresponding to space- like trajectories [18]. 

If such event distributions were verified by employing high-precision neutrino-speed 
experiments, the presented formalism would provide a direct means for experimentally 
determining the absolute masses of neutrino mass eigenstates. We add the following caveat 
to this conclusion: given that the canonical flavour oscillation formalism does not take into 
account the mass-dependent amplitudes for space-like separations, such a determination 
would require a generalisation of the standard neutrino oscillation phenomenology. 

3.3 Short baseline neutrino-speed experiments 

For short-baseline experiments, such as LSND [40, 41], the amplitudes (3.1) are only appli- 
cable to non-oscillation events and we will consequently restrict the discussion that follows 
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Figure 3. A graphical representation of the probabihty to detect a single mass eigenstate outside 
the light-cone. The ordinate aio corresponds to the ratio of the probability at T and the probability 
at Fio := 10 ns. 

to such events.^ For comparison with the above analysis, we consider m = 10~^ neV 
and find that a reduction of the basehne from ^op^ra _ 730 km/L^°^^^^ = 735 km to 
^LSND = 30in yields a suppression of ar(r). Interestingly, by setting m = 1 /xeV one 
finds that a reduction of the baseline to 30 m now results in a practically indistinguish- 
able ar(r). This can be understood by analysing the form of equation (3.1); the modified 
Bessel function of the second kind approaches zero for an increasing argument, so that the 
contributions to the integration of equation (3.1) over large values of r] become negligible. 
Therefore, for a sufficiently small mass the values over which the modified Bessel function of 
the second kind contributes to the integration is greater than the baseline of = 30 m. 

By increasing the mass, however, all the contributing terms to the integration can be shifted 
to within the considered baseline. 

4 SN1987a 

Repeating the above analysis for the scenario of SN1987a by adjusting the baseline to 
^sNi987a ^ ^ g ^ 10^ ly [42-44] yields effectively indistinguishable results from those shown 
in figures 2 and 3, for their respective mass scales. This may once again by understood 

®The LSND neutrino-oscillation data requires a mass-squared difference of the order of an eV^. It is 
unlikely that this particular mass-squared difference is associated with the fieV- neV range mass eigenstates. 
If the latter mass eigenstates do exist, these will be detected mostly in the non-oscillations events. 
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by considering the form of equation (3.1): over large values of r/, the contribution to the 
integration by the modified Bessel function of the second kind becomes negligible. The 
significance of this is that, for the considered mass scales, any early arrival time due to this 
effect would be of the order of nano seconds, even for baselines of astrophysical magnitude. 
This is consistent with the data recorded for SN1987a. 

5 Conclusion 

A preliminary quantum field theoretic investigation of the possible detection of a parti- 
cle across a space-like separated interval was undertaken. We have found that for suf- 
ficiently light mass eigenstates space-like amplitudes become significant for macroscopic 
length scales. Furthermore, we have shown that such amplitudes are only applicable to 
states with a vanishing momentum expectation value. Therefore, a formalism incorporating 
such amplitudes for a moving mass eigenstate was developed and applied to the OPERA 
and ICARUS neutrino-speed experiments. This yielded a mass dependent distribution with 
a maximal detection probability centred around a light-like trajectory. 

For comparison with the OPERA experiment a choice of m ~ 10~^ neV yielded a sig- 
nificant detection probability for events with an early arrival time of 60 ns; however, given 
that the predicted event distribution is maximal on the light-cone it did not reproduce 
the peak at 60 ns of the OPERA event distribution. In order to compare the developed 
formalism with the ICARUS experiment a choice of m ~ 1 ^ueV yielded a visual agree- 
ment between the predicted and the observed event distributions. Given that the currently 
available data consists of only seven data points it was, however, not possible to calcu- 
late a statistically meaningful fit. Nevertheless, this raises the possibility of employing 
high-precision neutrino-speed experiments to determine the absolute masses of neutrino 
mass eigenstates. Furthermore, it was shown that the predicted space-like amplitudes are 
consistent with the empirical data on SN1987a. 

A Notational details 

In equation (2.1) a takes on the values ±1/2, and m refers to the mass of the particles 
and antiparticles described by the field ^(x); is the time component of the energy- 
momentum four vector p'*, and equals y^p-^ + . Any such can be expressed in terms 
of the standard four-momentum := (m, 0) by the usual Lorentz transformation: 



The expansion coefficients given in equation (2.1) differ from those given in reference [10] 



(A.l) 



by a factor of \/rnJjP, such that 



u{p,a) = D{L{p))u{Q,a) 



(A.2) 



v{p,a) = D{L{p))v{Q,a) 



(A.3) 
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with 



D{L{p)) 



E + m 1 + 

I E+m 

2m \ O 



1 



crp 

E+m, 



(A.4) 



Here, cr := (ax, CTy, az), are the standard Pauh matrices; 1 and O are the 2x2 identity and 
nuh matrices, respectively. The expansion coefficients at rest are defined by 



«(0,l/2) 



^;(0,l/2) 



1 



1 



^(0,-1/2) 



v{0,-l/2) 



1 

71 



1 

71 



(A.5) 



(A.6) 



where the index i introduced in equation (2.1) is the row index of the above defined 
expansion coefficients. In reference to equation (2.3) 



1 



7 



1 



(A.7) 



Furthermore, the vacuum |vac) is Lorentz invariant with (vac | vac) = 1, and the annihilation 
and creation operators satisfy the usual fermionic anticommutation relations 



[a\p',a'),a\p,a)}=0 
{p',a'),a{p,a)^ =0. 



(A.^ 



In equations (2.2) and (2.4), * represents complex conjugation and T denotes matrix trans- 
position. The following important convention is then used 



a*(p,cr) = a^p,a). 

The counterpart of equations (A. 8) and (A. 9) holds for b{p,a) and b'^{p,a). 



(A.9) 
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